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2011 Æ 12 A 11 日 是 西安 交通 大 学 杰出 校友 钱学森 先生 
的 百年 诞 展 。 为 缅怀 钱学森 学 长 ,学习 他 的 科学 思想 和 卓越 风 
范 , 展 示 其 丰功伟绩 和 人 格 魅 力 ,西安 交通 大 学 举办 了 * 纪 念 钱 
SRNE 100 周年 系列 活动 : 作为 制 片 方 之 一 ,参与 西部 电影 
集团 摄制 传记 故事 片 (钱学森 )》; 与 中 央 电视 台 合作 ,出 品 纪录 片 
《实验 班 的 故事 一 一 沿 着 钱学森 走 过 的 路 ); 扩 建 钱学森 生平 业 
绩 展 馆 , 向 校内 外 开放 ;举办 钱学森 科学 与 教育 思想 研讨 会 ;出 
版 发 行 (钱学森 力学 手稿 )《 钱 学 森 年 谱 ( 初 编 )》《 钱 学 森 第 六 
次 产业 革命 思想 探 微 丛书 ) 等 。 

钱学森 先生 在 美国 深造 和 工作 期 间 留 下 大 量 珍贵 手稿 ,这 
些 手稿 真实 展示 了 钱学森 先生 博大 精深 的 学 识 、. 开 拓 求 实 的 精 
神 和 严谨 奋进 的 作风 ,是 钱 老 勇 攀 科 学 高 峰 和 严 谭 治 学 的 集中 
体现 。 这 里 ,我 们 将 部 分 原稿 整理 汇集 成 册 , 出 版 (钱学森 力学 
手稿 ), 作 为 钱 老 百年 诞 展 的 献礼 。 | 

《钱学森 力学 手稿 } 共 10 卷 ,包含 两 部 分 内 容 。 第 一 部 分 是 
草稿 ,包括 扁 壳 、 球 壳 和 圆柱 壳 届 曲 分 析 的 公式 推导 和 数值 演 
算 。 在 研究 圆柱 壳 轴 压 届 曲 问题 时 ,为 了 求 得 圆柱 壳 体 的 临界 
压力 ,在 有 关 的 五 百 多 页 草稿 中 ,对 多 达 二 十 多 种 可 能 的 届 曲 模 


态 逐 一 进行 公式 推演 和 数值 计算 ,最终 才 找到 满意 的 并 在 论文 
中 采用 的 届 曲 模 态 。 仔 细 观 察 草稿 中 的 数据 列表 ,每 个 数字 有 
效 位 数 都 长 达 八 位 ,在 手 播 机 械 式 计算 机 作为 主要 计算 工具 的 
年 代 , 这 串 串 数字 凝聚 着 多 少 现今 难以 想象 的 艰辛 劳动 。 

第 二 部 分 是 手稿 ,以 航空 航天 工程 为 核心 ,涵盖 空气 动力 
学 、 固 体力 学、 火箭 技术 .工程 控制 论 和 物理 力学 等 领域 的 部 分 
学 术 论文 手稿 .打印 稿 和 讲义 。 

《钱学森 力学 手稿 ) 是 在 西安 交通 大 学 校 领导 的 大 力 支持 
下 ,由 西安 交通 大 学 航天 航空 学 院 沈 亚 鹏 教授 整理 完成 。 图 书 
出 版 过 程 中 得 到 了 西安 交通 大 学 党 委 宣传 部 .校友 关系 发 展 部 、 
图 书馆 、 航 天 航空 学 院 等 的 积极 协助 ,在 此 深 表 感 谢 . 
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EMISSIVITY OF DIATOMIC GASES AT LOW PRESSURES 


I. INTRODUCTION 


30 July 1951 


Emissivity calculations for diatomic gases from spectroscopic data 


were developed recently by S. S. Penner (Ref, 1). His method is based 
upon the use of an average absorption coefficient for the entire funda- 


mental and higher vibration-rotation bands, The method is thus effective 


when there are extensive overlapping and broadening of the spectral 


lines, and hence is accurate for gases at high total pressures and 


temperatures. At low pressures, the lines do not overlap and a different 


approach to the problem should be made, Penner and M, H, Ostrander 
(Ref. 2) have computed the emissivity of carbon monoxide for the case 
of non-overlapping lines by a numerical procedure, using spectroscopic 
data obtained recently by Penner and D. Weber (Ref. 3). The results are 
in excellent agreement with the emissivity determined experimentally by 
W. Ullrich and H. C. Hottel (Ref, 4). The amount of numerical work 
involved is, however, rather heavy, It is the purpose of the present 
paper to develope an avproximate but convenient formula for calculafing 


the emissivity of diatomic gases for the case of non-overlapping lines. 


II. FORMULATION OF THE PROBLEM 


If T is the temperature, @ the characteristic temperature, 


Y the wave number, y* the characteristic wave nu:ber, P, the spectral N a 
JA coefficient at » , p the partial pressure of the emitting gas, 
4): UA)! 
7 A y YE 


figo? 


sonale 


and L the optical path length, then the emissivity £ of the gas under 
the specified conditions is 


fofo! fare 
/ n de: / KZ , (1) 


۶ * > A, dran fm 
«A Cum Longo olay Ma E AA fb ul 
If only the fundamental vibration-rotations band is idered, the mall 
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absorptions for the lines centering on the wave numbers corresponding 
to the indicated transitions. The Ss! can be computed in turn by using 
the results of J. R. Oppenheimer (Ref. 5). As 
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where N, is the mmber of emitting molecules at tempetature T per jum sort. | 
volume per unit pressure, © the effective enan, dl the reduced mass, wee 1 of à 
e the velocity of light, Q the complete internal partition funeti 

Ets are the internal energy levels given by LS 
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where x, J >£ are molecular constante in their standard notations, 
These constants are non-dimensional and are small, The F's and G's are 


FI), VE E i, 


Lo E -4 < TP Zë i 
Ey ۶,۷ bi li Cp L 


6 < | + w |- Ag 


4 1- TEE GC 


og 4 lue ` dea el ap^ 
. 
n ا‎ s ige 

4 


pm. 2 
q = Af. ali Hc d Gap 
If the fundamental vibration-rotation band gives the main con- 
tribution to the emissivity of the gas, the above equations give the 
necessary information to calculate approximately the emissivity ¢ . 
III. APPROXIMATE SOLUTION 


The numerical work is carrying out the computation indicated in the 


previous section is very heavy. A short formula, however, can be 

La 
developed: First of all, when the lines are from each other, 
each line can be considered alone, independent of others. Furthermore, 

4 

the value of the factor outside of the bracket in the mmerator of Eq. )32( 
can be approximated by its value at the center of each line, Thus 
according to S. S. Penner (Ref, 6) 
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where the Ps are the absorption coeffieient due to the particular line 
with transitions as indicated, The integrals can be easily evaluated 
(Ref. 6) and are given by the modified Bessel functions I, and D: 


+ (1- Ae AB) = raf) 5e JE (5) ZAA (10) 
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A further approximation can now be made, The magnitude of 
T's and ) is are generally quite large if the product pL of pressure 
and optical path length is of the order of unity, Therefore, the 
asymptotic values of the Bessel functions can be used, Then 
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By substituting Eqs. (14) and (15) into (9), the emissivity is calculated 
as a sum over j, . 

To carry out the sum over j, one can use the Euler-Maclawrin 
sumation formula (Ref, 7), which evaluates the sum by an integral, 
First, due to the smallness of , £ , the following expansions, 
including terms up to the square of A , are approximate, 
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The corresponding quantities for the transitions j - 1 -> j can be very 
easily obtained from Eqs. (16) to (20) by replacing j with -j. [heane 
of this property of symmetry, the sum of terms from the transition 

j — j - 1 and the transition j - 1 -> j for every j is a function of 


32 only, Thus, VV Q AA. tems, 
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A is thus a constant independent of temperature and pressure, The f 
function is simply deduced from the ion Q as given by 
Eq. (8): 
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from the expansions given in Eqs. (16) to (dè. It y e 
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the sum in the emissivity E . The resulting integral over j extends 
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simply deducting the approximate value of the integral from 0 to 1 from 
the extended integral, Thus 
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The PTG) has the numerical value of 1,225, Finally then the expression 
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where f and g are functions given previously in Eqs, (23) and (24). 
Since the value of / is nearly unity and the factor before 
g in Eq. (9 is small, a good approximate equation for t he emissivity 


is 
> “- Y 
EN TRF 
IV. APPLICATION TO CARBON MONOXIDE 


For carbon monoxide, the molecular constants are 
b= 6,۲ 
y'- 2142.3 a 
Y 895 
$ * 0,0091 
X = 0 
The value of A computed from the measurements of Penner and Weber 


(Ref. 3) 1s 22.95 
3 
A = Dafra! e" 
They have also determined (Ref, 8) b to be 0.077 ca at one atmosphere 
of total pressure, 


10011111110 
T = 300°K and a total pressure of one atmosphere 


۵2+ 10 FV pL 
1.593 
E ms x 10? ۳ pL 
where pL is in atm-cm. By using the more exact equation (25), the 
enissivity is i 


t 
£ = هید‎ x10? Jet 
Va? l 
The difference between the approximate value and the more exact value is 


quite small, The comparison between the computed emissivity and the 


measurements of Ullrich and Bette) (Ref. 4) is shown, in Fig. 1, The 


agreement is quite satisfactory up to pL of approximately 10, 


as 


2. 


3. 


4. 


5. 
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Summary 
This paper shows that the combustion in the rocket motor can be stabilized . 
against any value of time lag in combustion by a feedback servo link from a 
chamber pressure pickup, through an appropriately designed amplifier and é fe a 
| control capacitance, aetémw on the propellant feed ine. The technique of 
stability analysis is based upon a combination of Shono diagram and Nee 
diagram. Jw iii o d del m, irá Las fue 3 
nafa dallo t solih "en me eee Mer, He ¿saa $ 
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* Robert H. Goddard Professor, Deniel and Florence Guggenheim Jet Propulsion 
Center, California Institute of Technology. È 


The phenomenon of rough burning in liquid propellant rocket motor has 
been interpreted as the instability of the coupled system of propellant feed 


and combustion chamber by D. F. Gunder and D. R. Friant (Ref. 1), M. Yachter 


(Ref. 2), .لا‎ Summerfield (Ref. 3), and L. Crocco (Ref. 4). The essential 
feature of these theories is the time lag between the instant of injection of 
the propellant and the instant when the propellant is burned into hot gas. 
Crocco has further improved op this concept by considering the time lag as an 
integrated effect of consecutive stages, each of which is controlled by tho 
prevailing pressure in the combustion chamber. As a result of thia new concept, 
Grocco showed the possibility of intrinsic instability with constant injection 


rate not influenced by the chamber pressure. 

The present paper will first give a slightly more general formulation of 
Oroeco's concept of time lag, allowing arbitrary pressure dependence of lag. 
Then the problem of intrinsic stability is discussed by applying e method 
suggested by M. Satche (Ref. 5). This method is based upon a modification of 
the Nyquist diagram and is particularly useful for systems having time lag. For 
easy reference, this now diagram will be called the Satche diagram. The later 
sections of the paper will show the possibility of stabilizing the combustion 
by means of a feedback servo for all values of time lag. Such possibility of 
servo=stabilisation was first mentioned by V. Bollay in his admirable paper 
(Ref. 6) on the application of servo-mechanisus to aeronautics. The present 
study definitely shows the power of this idea. 


Time Lag in Combustion 

Let m(t) be the mass rate of generation of hot gas by combustion at time 
instant t. Consider, for simplicity, a monopropellant motor. Then the mass rate 
of injection at t can be denoted by mq (t). Let (t) be the time lag for that 


paréel of propellant which is burned at the instant t. Then the mass burned 
during the interval from t to t + dt must be equal to the mass injected during 
the time from t=- to te + d(te ). Thus 


a) 


The mass of hot gas generated is either used to fill the combustion chamber 
by raising its pressure p(t) or is discharged through the rocket nozzle. If the 
frequency of the possible oscillations in the chamber is small, then the pressure 
in the chamber can be considered as uniform and as a first approximation (Ref. 7) 


the rate of flow through the nozzle can be taken as proportional to the 


instantaneous chamber pressure p(t). Thus if is the steady mass rate flow 
through the system, is the average mass of hot gas in the chamber, and if 
the volume occupied by the unburned liquid propellant is neglected, 


(2) 


where 
t y gon ñ By following Crocco, the non-dimensional variables for the chamber pressure 
KE and the rate of injection are defined as 


is the steady state pressure in the combustion chamber. 


(3) 


Ze: and are then the fractional deviation of pressure and injection rate 
i from the average. With Eq. (3), can be eliminated from Eqs. (1) and (2), 
7 Ñ and 


(4) 


! fo calculate the quantity » Crocco's concept of pressure dependence 
K » 14) ; of time lag has to be introduced. If the rate at which the liquid propellant is 


prepared for the final rapid transformation into hot gas is a function 
then the lag is determined by 


= Constant 


By differentiating Eq. (5) with respect to t, 


The concept of small perturbation from the steady state will now be explicitly 
introduced: assumo that the deviation of the pressure p from the steady state 
value is small. Then at the instant t and at the instant 

te can be expanded as Taylor's series around +. By taking only the first 


order terms, 


Here 1s the lag at the average pressure , a constant now. Then 


By combining Eqs. (4) and (6), the following equation is obtained 


Sa 
and 
(9) 
If n is a constant independent of » then is proportional to ° 
This is the form of assumed by Crocco. The present formulation of the 
problem is slightly more general in that is arbitrary and the value of 2 
is to be computed by using Eq. (8), and is a function of a is, of course, 


the gas transit time. 


Intrinsic Instability 
Crooco called the instability of combustion with constant rate of injection 
the intrinsic instability. If the injection rate is constant not influenced by 


the chamber pressure p, then e Therefore the stability problem is controlled 


by the following #imple equation obtained from Eq. (7), 


(10) 


Now let 


Then (11) 


This is the equation for the exponent s. 

Crocco determined the value of the complex number s by studying the set of 
two equations for the real and the imaginary parts of Eq. (11). However if the 
point of interest is whether the system is stable or not, one can use the well- 
known Cauchy's theorem with advantage. Let 


(12) 


Then the question of stability is determined by whether has zeros in 
the right half of the complex s-plane. This question itself can be in turn 
answered by watching the argument of when 8 trace a contour enclosing 
the right half seplene. Specifically, let s trace clockwisely the contour 
consisted of the imaginary axis and a large half circle to the right of the 
imaginary axis (Fig. 1). If the vector M a number of complete 


clockwise revolutions, then that number is, according to Cauchy's theorem, the 


difference between the number of zeros and the number of poles of in the 
right half seplane. Since evidently has no poles in the s-plane, the 
number of revolutions of is the number of zeros. Hence for stability, 

the vector must not make eny complete revolutions, as s traces the 
specified contour. Therefore the stability question can be answered by plotting 
graphically on the complex plane. This graph is, of course, the well» 
known Nyquist diagram. 

A direct application of this method to given by Eq. (12) is however 
inconvenient for the complication caused by lag term (Ref. 8). Me Satche 
(Ref. 5) however proposed a very elegant and ingenious method to treat such 
system with time lag: Instead of „ break it into two parts, 


(13) 


The vector is thus a vector with vertex in and its tail on 5 
The graph of is the unit circle for s on the imaginary exis. For a on the 
large half circle, is within the unit circle. The graph of is 


Im 


the straight line (Fig. 2) paralled to the imaginary axis when s is on the 
imaginary axis. When s is on the large half circle, is a half of 
a great circle closing the contour on the left. A moments! reflection will 
show that in order for the vector not to make complete revolutions 
for any value of „ the contour must lie completely out of the 
contour. That is, for unconditional intrinsic stability, 


or (15) 


` 


When , the contour and the contour intersect. Stability 
is still possible however, if for within the unit circle (Fig. 3), is 
to the right of e This condition is satisfied if 


Or if 


where (16) 


When » then with 


Q7) 
« Therefore when " has the oscillatory solution with the 
angular frequency  *. 

These results on intrinsic stability were obtained by Crocco. The present 
discussion with the Satche diagram however seems to be simpler. For the more 
complicated stability problem treated below with feed-system and servo-control, 
the solution is hardly practical without the Satche diagram. 


m 3 El 
Consider now a system including the propellant feed and a servo-control 


represented by Fig. 4. In order to approximate the elasticity of the feed line, 


Y 
y 
yo 
li 
* 
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a spring lead capacitance is put at the midway point between the propellant 
pump and the injector. The spring constant is to be computed from the feed 
line dimensions." Near the injector there is another capacitance controlled 

by the servo. The servo receives its signal from the chamber pressure 

pickup through an amplifier. If the feed system and the motor design is fixed 
by the designer, the question is whether it is possible to design an appropriate 
amplifier so that the whole system will be stable. Because there is no accurate 
information on the time lag of combustion, a practical design should specify 
unconditional stability, i.e., stability for any value of . 

Let be the instantaneous mass flow rate out of the propellant pump 
and be the instantaneous pressure at the outlet of pump. The average flow 
rate must be e The average pressure is e The pump characteristics can 
de represented by the following equation, 


(18) 


If the time rate of change of mass flow is small, is simply related to the 
slope of the head-volume curve of the pump at constant speed near the steadye 
state operating point. For constant pressure pump or the simple pressure feed, 
is zero. For conventional centrifugal pumps, is approximately 1. For 
displacement pumps, is very lerge. 
Let be the instantaneous mass rate of flow after the spring loaded 
capacitance, be the spring constant of the capacitance and the 


instantaneous pressure at the capacitance. Them 


(19) 


* See the Appendix for details. 
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where is the deneity of the propellant, a constant. 

In the following calculation, the pressure drop in the line by friction:l 
forces will be neglected. Then the pressure difference is due to the 
acceleration of the flow only. That is 


(0) 


where A is the cross-sectional area of the feed line, a constant, and is the 
total length of the feed line. Similarly, if is the instantaneous pressure 
at the control capacitance, 


If the mass capacity of the control capacitance is C, then 
(22) 


Since the control capacitance is very close to the injector, the inertia of the 
mass of propellant between the control capacitance and the injector is 
negligible. Then 


(23) 


where A, the effective orifice area of the injector. Au can be eliminated 
from the calculation by noting that at steady state, the difference of pressures 
and ‚or is 


(24) 


Equations (18) to (24) describe the dynamics of the feed system. By a 


straightforward process of elimination of variables, a relation between 


end C is obtained. To express this relation in non-dimensional form, the 
following quantities are introduced, following the notation of Crocco: 


(25) 


and (26) 


where the gas transit time given by Za. (9). Then the non-dimensional 
equation relating 5 and is 


where is the non-dimensional time variable defined by Eq. (9). 

The dynamics of the servo-control is specified by the composite of the 
instrament characteristics of the pressure pickup, the response of the amplifier 
and the prepértios of the servo. Since it is not the purpose of the present 
paper to discuss the detailed design of the servo-control, the overall dynamics 
of the sorvo-control is represented by the following operator equations 


(28) 


whore Y is the ratio of two polynomials with the denominator of higher order 
than the numerator. 


Equations (7), (27) and (28) are the three equations for the three 
variables 5 and . Since they are equations with constant coefficients, 
the appropriate forms for the variables are 


W ai 
ME 


ati 


(29) 


By substituting šq, (29) into Eqs. (7), (27) and (28), three homogeneous 
equations for a, b and c are obtained. In order for a, b, c to be non-zero, 
the determinant formed by their coefficients must vanish. This condition can 


be written as follows: 


(30) 


This is the equation for determining the exponent s. F(s) is now recognized as 
the overall transfer function of the servo-control link. The complete system 


stability depends upon whether Eq. (3) gives roots that have positive real part. 


Y 5 
The system characteristics without the servo-control can be simply obtained 

from the basic equation (30) by setting F(s) = 0. Let it be assumed that the 

polynomial multiplied into has no zero in the positive half s-plane, as is 

usually the case. Then Eq. (30) can be divided by that polynominal without 

introducing poles in the positive half s-plane into the resultant function. 

That is, for the Satche diagram, one has again 


is thus again tho "unit ۰ 1s now much more complicated: 


-12- 


(31) 
The intercevt of when s is pure imaginary is given by setting s = 0 in 
GP (31), 1.000 

(32) 
Since all the parameters ` „ P are positive, the magnitude of is 
now smaller than the magnitude of given by Eq. (14) for the intrinsic 
stability problem. Thus the effect of the feed=system is to move the 
curve towards the unit circle of in the Satche diagram. For instance, for 

" is just tangent to the unit circle for the intrinsic system 


without considering the propellant feed. But with the propellant feodesystem, 
contour will intersect the unit circle and the system will become unstable 
for time lag exceed a certain finite value. The influence of the feedesystem 
is thus always de-stabilizing. This is further confirmed by considering the 
asymptote of for large imaginary » obtained from Eq. (31). That is 


(33) 


Therefore for large imaginary 8 , approaches asymptotically a line 
parallel to the imaginary axis at a distance 


to the left of the imaginary axis. The effect of feed-system is again to nove 
towards the unit circle. 


1370 
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It is thus evident that for the parameter near 1/2 or larger than 1/2, 
it would be impossible to design the system for unconditional stability. In 


the Satche diagram, contour and will always intersect without a 
servo-control. 
Complete Stability with Servo-Control 


If the polynomial H(s), 


(34) 


which multiplies into in Eq, (30), has no poles and zeros in the right 
half e-plane, then the occurance zeros of the expression in Eq. (30) in the 
right half s-plane cen be determined from the Satche diagram with 


and 

(35) 
As a traces the contour of Fig. 1, is again a unit circle. Therefore 
1f simultaneously the contour is completely outside the unit circle, 


there can be no root of Eq. (30) in the right half s-plane, In other words, if 
the transfer function F(s) of the servo-control link is so designed as to place 
the contour completely out of the unit circle (Fig. 5), then the system 


is stabilized for all time lags. 


Then without the servo-control, the 


Of primary interest is the behavior of when 8 is a pure imaginary number 
05 real. Thus 


This contour for is plotted in Fig. 6. It is evident that for sufficiently 
large values of time lag, the system will be unstable. On the other hand, if 
the contour can be changed by the servo-control to say 


Then as plotted in Fig. 6, the new contour is completely outside of the 
unit circle of « Therefore the system is now unconditionally stable. A 
straightforward calculation from Eqs. (31) and (35), shows that the required 
transfer function F(s) for the servo link is 


The servo link has thus the character of an integrating circuit. If with given 


response of the chamber pressure pickup and of the servo for the control 


capacitance, en amplifier could be designed to give an overall transfer function 
close to that specified above, the combustion cen be stabilized by such a servos 
control. 

As the second example, take 


Since » the feed pressure is thus constant with even variable flow 
of propellant. The case then corresponds to that of a simple pressure feed. 
Without the servo-control, 


When s is pure imaginary, 


This contour of is plotted in Fig. 7. It is evident that without servos 


control the combustion will de unstable for sufficiently long time lag. In fact, 


the system is even less stable than the system considered in the first example: 
It will become unstable at shorter time lag. The part of the contour near 
2 is of special interest. Near » the contour comes so close of 
the unit circle of that if the value of time lag is such as to make 
and for 2 very close to each other, then an almost undamped oscillation 
at 2 can occur. This critical value of 18 evidently smaller than 
the critical determined from the true intersection of with the unit 


-16- 


circle at 0.65. Such near instability at smaller values of time lag 
can be easily overlooked in the analytic treatment of the stability condition 
by Crocco, end yet such possible instability should not be dismissed. This, 
perhaps, indicates the superiority of the present graphical method. 

For unconditional stability, should be displaced out of the unit 
circle, to, say, the same "stable" contour as in the first example. The required 


transfer function F(s) is calculated to be 


The required servo link must then have the character of double integrating 
circuit. Furthermore, the transfer function has two purely imaginary poles at 
» This unrezlistic requirement on the amplifier comes from the original 
feed-system dynamics and is due to the neglection of frictional damping in the 
feed line. In any actual system, the frictional damping in the feed line will 
remove these purely imaginary poles of the required transfer function F(s) and 


replace them by two complex conjugate poles. 


Stability Criteria 
In the preceding discussion of servo-stabilization, it is assumed that the 


polynomiel H(s), Eq. (34) has no pole or zero in the right half s-plane. This 

is however not necessarily the case. In general then, one should first investigate 
the number of zeros and poles of H(s) in the right half s-plane. To do this, it 
should be recognized that the polynomial in Eq. (34) before the factor F(s) 
usually does not have — in the right half s-plane. Therefore instead of 


studying H(s), one can study the ratio of H(s) and that polynomial. That is, 


2 ace OAN” 2 — a 
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the number of zeros and poles of H(s) in the right half s=plane is the sane 


es the number of zeros and poles of the following function 


(36) 


where 


(37) 


According to the Nyquist criterion, the number of poles and zeros for 1+K(s) 

in the right half s-plane can be found by plotting the Nyquist diagram of 1+K(s) 
with s tracing the contour of Fig. 1. In fact, 1f 1+K(s) or H(s) has zeros 
and poles in right half s-plane, then K(s) will carry out clockwise 
revolutions around the point =l, as s traces the contour of Fig. l. Hence the 
necessary information on H(s) can be obtained by plotting the Nyquist diagram of 
KLS), 

When one divides the Eq. (30) by H(s) in order to obtain and as 
given by Eq. (35), zeros and poles ere introduced in the right half 
seplane. The poles of K(s) mst come from F(s), since the polynomial in the 
denominator of Eq. (37) has no zero in the right half seplane. Therefore the 
original expression in Xq. (30) also has poles in the right half s-plane. 
Hence in order for the original expression in Eq» (30) to have no zero in the 
right half seplane, must make clockwise revolutions 
around the unit circle, In order the stability be unconditional, i.e., stable 
for all time lag, the contour should never intersect the unit circle. 


Therefore the general unconditional stability criteria are, first, contour 
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completely outside of the unit circle; and, second, making counter 
clockwise revolutions around the unit circle as s traces the conventional 
contour enclosing the right half s-plane, These are the criteria for stability 
with the Satche diagram. To determine » one has to use the Nyquist diagram 
of K(s), Eq. (27). Thus the stability problem for the general cose requires 
both the Satche diagrem and the Nyquist diagram. (Fig. 8) 


Soncludins Remarks 
In the previous sections of this paper, the theoretical possibility of 


completely stabilizing the combustion for any value of time lag by servo=control 
is demonstrated. The great flexibility of electronic amplifier seems to 
indicate that this theoretical possibility cen be always realized. On the other 
hand, without the servo link, unconditional stability is shown to be generally 
impossible. Therefore the concept of feedback servo is indeed a powerful tool 
in controlling the behavior of a timeelag system. It is to be realized, of 
course, that the proposed scheme is but one among many. No attempt is made here 
to give an exhaust treatment of all possible schemes. The best scheme is 
certainly to be determined by detailed considerations on all aspects of the 
engineering problem, such as the possibility of high frequency accoustic 
oscillations which are not considered here. The main purpose here is to give a 
general discussion of the concept together with a suggested general method of 
analysing the stability by the Satche diagram. 

It is of interest to point out that stabilization by servo=control is only 
one phase of the general concept of feedback link. The opposite case of 
de-stabilisation could be of importance also. For instance, consider the 
so-called valveless pulsejet. It is not always possible tc operate the engine 


with the desired pulsation. With a feedback servo linking the combustion chamber 


i 1 7 m 
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pressure pickup through an amplifier to the fuel line, the system can be dos 


stebilized at the desired operating frequency end thus operete the engine at 


that frequency of pulsation. This application of servo=destabilization gives 
the velveless pulsejet e new flexibility and an extended range of operation. 
Therefore 1% ۵ worthwhile to explore carefully all possible applications 
of feedback control to systems with time lag. 


Appendix 
a of e 


If L* and c* are the characteristic length and the characteristic velocity 
of the motor, and if To is the chamber temperature, R the gas constant, the 


transit time is 


To calculate J and E defined by Eq. (25), it is more convenient to use the 


average propellant velocity in the feed line. Thus 


Thus, according to Eq. (25) 


A consistant set of units would be in slugs per cubic foot, in feet 
per second, in feet, in seconds and in pounds ver square foot. 
If à the diameter of the feed line, h its thickness and Zt the Youngs 
modulus of the tube material, then s the chenge in volume of the feed line 


per unit rise in pressure, is 


Therefore Eq. (25) gives 


A consistant set of units would be in pounds per square inch; E?’ in 


pounds per square inch, in feet, in seconds, and in feet per second. 
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Contour traced by the variable s for the Satche Diagram or the 
Nyquist Diagram. 


Stable Satche Diagram for 1 

Unstable Satche Diagram for Intrinsic Oscillations; 
Servo-controlled Liquid Monopropellant Rocket Motor. 

Satche Diagram for the Original end for the Servo-Stabilized System. 


Satche Diagram for the Original end for the Servo-Stabilised System 
NS 


without servo intersects the unit circle; with servo 
is outside the unit circle. Numbers beside points are the value of 


Satche Diagram for the Original and for the Servo-Stabilized System 


P= 3/2, J= 4, E = 1/4, = 0, 


without servo intersects the unit circle; with servo 
is outside the unit circle. Numbers beside points are the value of 


Full curve for positive dotted curve for negative ۰ 
a) we Diagram for x(s), with two zeros for 1+X(S) in right 


half s~plane 
b) Corresponding 882516 Satche Diagram 
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Engineering Cybernetics 


1.19 Linear Systems of Constant Coefficients 
š Let us consider the simplest systern - a first order system. 
That is, the differential equation of the system is a first order linear 
differential equation of constant coefficients, If the system is assumed 
to be free and is not subjected to "forcing function, Mthen the differential 
equation can be written as 


É hy = 0 64. 1) 


4 may be called the spring constant and is real. When there is no 
variation of & with respect to time, 2 vanishes and then Eq. (1.1) 
requires Je O. Therefore the stationary state or the equilibrium state 
of the system corresponda to quo 3 

The solution of Eq. (1.1) is 


21% 


where A is the initial value of H or 


dá is thus the initial disturbance of the system from the equilibrium 
state. The behavior of the system for É>O , is illustrated in Fig. 1.1 
for both positive 4 and for negative 4 . It is seen that for £>0 ; 
the magnitude of decreases with time. Then as time increases 
indefinitely, onn . Therefore 


Fig. 1.1 


A 
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for ۳ م2‎ , the disturbance of the system will eventually disappear. 
The system can then be said to be stable. When ((o, the disturbed 
motion of the system increases with time and eventually the disturb- 
ance will become very large no matter how small the initial displace- 
ment is, and will never return to the equilibrium state once disturbed. 
Such systems are thus unstable. 

For systems of higher order, the differential equation will have 
higher derivatives. The n-th order system has the differential equation 


B ۳ 
+ ig =o an 


For a physical system, the coefficients u, و‎ ax are real. Then 
the solution of Eq. (1. 4) can be written as 


y= D P SF sin (aye + fa) (1.5) 


where % ٠ fh are real and are related to the coefficients By a da 
and A are the phase angles. The motion of the system is thus only 
stable if all Sa? are negative. If one of them is positive, the disturb- 
ance will eventually diverge, and the system is thus unstable. 

From the above examples it is seen that the crucial question to 
ask about the behavior of a linear system of constant coefficients is the 
question of stability. Needless to say, the usual aim of an engineering 
design is stability. The question of stability can be answered however 
once the coefficients of the differential equation are specified. In case 
of the simple first order system specified by Eq. (1.1), the only infor- 
mation that matters is the sign of the coefficient £ . 


ko Linear System with Variable Coefficients 
If there is a variable parameter in the system under study, the 


stationary or the equilibrium state of the system can be changed by 
changing this parameter. It is سيت‎ men to expect the coefficients of 
the linear differential equation describing the system to be also functions 
of this parameter. For instance, the aerodynamic forces acting on an 
aircraft are functions of the speed of the aircraft. H the speed of the 
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Chapter 11 
Method of Laplace Transform 


For linear differential equations with constant coefficients and 
with time É as the independent variable, the method of Laplace trans- 
form is particularly useful in finding the solution. Of course, the 
problem can be solved by a number of other methods; but the method of 
Laplace transform appeals specially to the engineering scientists in that 
it reduces all problems to a uniform basis. The procedure of solution is 
then standardized and a general approach is possible. The theory and 
practice of Laplace transform is discussed in many texts.* It is not the 
purpose of the present chapter to do this. The purpose here is rather to 
give a summary of results which are useful to our discussion in the sub- 
sequent chapters for easy reference. For details and proofs, the reader 
should consult the texts cited. 


2.1 Laplace Transform and Inversion Formula 
1£ get) is a function of time variable É defined for C .. then 
the Laplace transform Y(S) of ye) is defined as** 


so 
ves) = 4 e ge dt (2.1) 
o 
where S is a complex variable having a positive real part, £ »o 
For other values of S , the function Y(S) is defined by the analytic 79 


continuation. The dimension of Y(S) is the dimension of Y multi- 
plied by time. 

When Y(S) is known, the original function for which (5) is the 
Laplace transform can be obtained in all cases by the Inversion Formula: 


* See for instance, H. S. Carslaw and J. C. Jaeger. "Operational R 
Methods in Applied Mathematics", Oxford,(1941); or R. V. Churchill, 
"Modern Operational Methods in Engineering". McGraw Hill. (1944). | 


For morecomplete theory, one should consult G. Doetsch, "Theorie Eco 
und Anwe der Laplace-Transformation", J. Springer, Berlin. pout vè 
(1937); or 2. V. Widder, "The Laplace Transform", Princeton, ااا‎ | ST Ñ 


n We shall use throughout capital alphabet to denote the Laplace trans- ۳ 9 
form of quantities denoted by a lower case alphabet. : | 


3-5 
Therefore the error signal vanishes as - 
Consider now another example of the input: Let the input be 
sinusoidal, i 
£ e 
x(t) = 2 عفر‎ 
where X4 is the amplitude and w is the frequency. Then 
Me 
S) = — 
X( ) EG (3. 12) 


The output due to the initial condition is the same as before. The output 
due to input is given by 


/ 
isla Spam s A I O / 
Yo Kä (s-i) TS +1) | ک‎ + T zx | 
Therefore according to our dictionary, the output He is 
tv 4 5 - tr 
R ci at e 
yo) = - Pe ES EAS E 
/ 1 + Lwa, / + LWE, 


The first term is a pure subsidence and the second term is the steady 
state output. Thus me WE) 


RS 0 belə. 
= / 0 
a < Ju == FG) 


This is in full agreement with our general result given in Eq. (2.16). 
Since 


/ = / EA f 


—— (3. 13) 
/ + LWE, Vi + oP GIN SÉ, © d 


the steady state output can be expressed as 


¿[ut - fa (ou) J 


I se 
[ yt ld es È L 
“steady i+ 
Therefore the amplitude of the steady state output is reduced by the a 


factor Ay FE ENN in comparison with the input, and the phase of 
the output lags behind the input by the amount tan ( wè). For low 


ES WI TS LEO S 
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the aileron deflection Š . A" SR the roll angle ( is thus 


p= NIA s E 


peed, then the above equation becomes 
‘ee — 


. If the roll speed is zero at #=0 , then the transformed equation is 
(1s + 44) 9 - LA Foe 


The transfer function cs) is thus 


ا اھ ر E‏ 
= / هه | 4+ ACs) Is‏ 

The behavior of the system is determined by the transfer function is 
thus similar to the cantilever spring with dashpot and the simple lag 
network. Here the characteristic time @ is J / . If the damping 


is very small, then Z and the behavior of the system becomes that 
of the simple integrator. 


(3. 36) 


3.4, Second Order Systems 

p Let us return to the cantilever spring with a dashpot, Fig. 3.1. 
Only now we attach a mass n to the dashpot end, The mass will 
introduce an inertia force 22 20 , and the equation of motion is now 


m PER +47 a a te 


with the initial conditions 


40 =g 
e +- š | (3. 37) 


The differential equation of motion can be rewritten in more convenient 
form by introducing the following parameters: 


with G(s) given. The method of Evans determines such roots as 
functions of the gain 4 , and is thus called the root-locus method. 
When this is done, any set of specifications on the roots gives a proper 
choice of the magnitude of X . This method then goes much beyond 
the mere satisfaction of criterion a) of Section$2, but actually solves 
the design problem for all three criteria stated in that section. 

Now let G(s) be specified by its zeros 3 °° and 
its poles q , 1 . Then from the definition of gain given 


by Eqs. (3. 16), (3.21) and (3. 23), 


„5 — 
GG) = A eee, (4. 16) 


Ae CDE 
(-h.)(- fa) 
For physical systems, the polynomials in the numerator and the 
denominator of G(s) has real coefficients. Then the ps are either 
real or form complex conjugate pairs. Similarly the Ss are either 


real or form complex conjugate pairs. P 
y APT Á E PEE ۲ L Re 


. „J al a A uk 


ser = mal jo alo foni Ce. Ma thu, we thal consid A o 
de معدم"‎ nd faikee. وام‎ e d GI) a f quel نه‎ 
s tf uio Md ماد‎ alm, ne, EM. KA < supom + 
سس«‎ Pou 


he vector P Be 
7 goes from tos. 4 
^ E The vector Q. e 9. goes 


£ S 5 m, - 
om the to is the variable int in the com lex s ne 
By using Eqs. (4. 18) and ( 4. 19), G (s) can be wr itten as 


GOD = a dE 2. S e. 
(g ο يه)‎ e“ *) (Ge) 


ive GE mihe m 


— Hence we can write Eq. 


¿0 
Ces) Re 


R = A(RR--=- 2,/ € €, - - -- Rn) ا‎ 


and d 6 (gef. - rar mm 


Since P’s and Q's are magnitudes of vectors defined by Eqs. (4. 18) 
and (4.19), they are positive. Therefore 2 is positive. The basic 
equation for the roots of inverse system transfer function, Eq. (4.15), 
is thus -40 

* e =-/ 

KR 
Therefore to satisfy this equation, we must have 


KR =! 


and 6 = 


The method of Evans consists of two steps: The first step is to ۰ 
determine all S that satisfy the appropriate angle condition of Eq. (4. 24). 
Then knowing such root-locus, we can compute Æ and hence £ , by 

Eq. (4.23), for each point on the root-locus. Evans has developed a 
number of useful rules for plotting the root-locus. We shall explain 

these rules presently. 


Rule 1 For Á =0, Eq. (4.15) shows that G(S)>oo. Thus for KO. 
the roots of VES are poles of G(s) , or the root-locus starts at the 


oscillating servomechanisms are less flexible than are oscillating 
control servomechnisms in which the oscillation is supplied by an 
independent generator. 
An elementary precaution to be observed, in order that the 

curve, which is constrained to pass through the point -2 , shall -void 

the neighborhood of the point -/, is that the curve shoul.{ interest the 
real axis, at the point -2 , perpendicularly. This imvlies thatthe vector 

/ / (éd) should be varying slowly in amplitude, and rapidly in angle, at 

the frequency at which the system oscillates. 


A relay . is a nonflinear device. But by mixing the 
signal with a sinusoidal oscillation gf high frequency and large amplitude, PA * 
the output is made to be linear with |respoct to the signal. Thus the T ۱ 
essential concept of oscillating control se chanisms is the lineariz- 
ation of a non-linear system. J. M. * has shown that this concept is 
applicable to any non-linear system, anc he calls this method the general 
linearizing process for non-linear contro. systems. We shall call the MS 
resulting servomechanism the general oscillating control servomechanism. 

Let us consider a general function 4%) where £ i the output H 
and = is the input. If instead of the variible x , we substitute the 
sum X+€ where € is much smaller than x . Then if the function Jeu 
is regular, we can expand jere into a Taylor series as 


2 
fore ge E) ék) 5 (6. 15) 


We now specify the input x as a periodic function of time É with the Ze 

period 7 , and E asa constant. Thef it دز‎ clear that Je is also 1 

a periodic function of time with the same period 7 . Some is true for J 
fd and d Y/ de . Periodid furctions can be expanded into 

Fourier series; thus if Së neglect powers uf € higher than first, we have 


pre Z 4, +> وی‎ +h, Sin not) 
rela 5 (a, Cos nut + b Sin nut )| (6. 16) 


ULTI 


* J. M. Loeb, Annfiles des Télécommunications, 5:65-71 (1950). 


We proceed to evaluate the right-hand member of (7. 11) by the method 


of residues. 

The integrand has certain poles, the poles of Z (s) . lying to the 
left of the path of integration, and other poles, which are the roots of the 
equation Sé SS 20, lying to the right of the path of D. npe m 
is easily seen that the integration upward along the line £ - too to js L 
is equivalent to integration in the clockwise direction along the closed 
contour formed by that line and the infinite semicircle in the right half- 
plane. Hence the right-hand member of Eq. (7.11) is -É times the 
sum of the residues of the integrand with respect to the several roots of (Ek 
the equation /- REL a el 

Now the typical root of the equation is ge $ + RIM go where m E 
is an integer, and the residue of the intagrand with respect to that pole 


is.— z£ (s+ en 8 Therefore finally 


E fs) E: (sani) 1 u 1.10) 


= ~o 


+ 
This formula ۳۹ considerable insight into the properties of E 6) . 
and at times may be useful in making approximate caleqigtions How- 
ever, we can easily obtain an exact representation of bet (s) in finite 


form. 
The function E (s) can be represented as the sum of a finite 


number of partial fractions, thus: 


». 
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Chapter VIII 


Linear Systems with Time Lag 


In this chapter. we shall introduce another new element into our 
linear systems with constant coefficient: the time lag. By time lag, سج‎ 
we mean that the relation between the different variables of the system 
cannot be expressed as a relation of these variables all taken at some 
time instant ¢ ; but on the contrary, the relation involves variables. 
some taken at the time instant É , and some taken at an earlier instant 
t-T. Those taken at the instant E- F then lag by the interval C behind 
the variables taken at the instant É . This time lag is thus quite different 
from the characteristic time constant of a first order linear system 
introduced in Section 3.1. Time lag systems are represented by 
differential-difference equations of constant coefficients and are more 
complex than the linear systems studied previously which are represented 
by differential equations, Systems with time lag were studied by many 
investigators; for instance, A. Callander, D. Hartree, and A. Porter* 
and N. Minorsky.** Our interest here is, however, somewhat more 
restricted. We wish to know: How can we analyse the performance of a 
feedback servomechanism if there is a characteristic time lag 2 inthe 
system? We wish, epecifically, to modify the method of Nyquist of 
Section 4. 3, to time lag systems. 

We shall develop§ the theory by treating a particular example of 
such systems; the example of stabilizing the combustion in a rocket motor 
by feedback control. The problem of combustion instability in rocket 
motors was treated by many authors, the following analysis of combustion 
lag time originates from the work of L. Crocco. ### For simplicity of 
calculation, **** we shall consider only the case of so-called low frequency 
oscillation in a rocket motor using single liquid propellant. 


* A. Callander, D. Hartree and A. Porter, Phil. Trans. Royal Society 
of London (A), 235:415-444 (1935). 


** N. Minorsky, J. Appl. Mechanics (ASME) 9:67-71 (1942). 
ses L. Crocco, J. American Rocket Society, 21: 163-178 (1951). 


d The following discussion is based upon a paper in J. American Rocket 
Society, 22:256-262 (1952). 
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Chapter IX 
Linear Systems with Stationary Random Inputs 


In the previous chapters, the inputs to a system are considered ۱ 
to be definitely specified functions of time É . However, there are 1 
many engineering problems for linear systems with constant coefficients 
where the inputs cannot be so definitely described. An example of such 
engineering problem is the problem of the motion and the stresses 
induced in the structure of an airplane wing in turbulent air stream. 
Here the input can be considered to be the time varying air-flow pattern. 
But the airflow pattern cannot be described as a definite function of time, 
but has to be recognized as a random function of time, specified by 
certain statistical characteristics. It is then evident, the output of the 
system, the stresses in this case, must be also a random function and 
can be described also only in statistical terms. The first objective of 
this chapter is then to find a convenient method of calculating the statis- 
tical properties of the output from the specified statistical properties of 
the input. This forms an easy extention of the early investigations by 
P. Langevin of the Brownian motion. 
Another example of random input is the so-called noise in control ited. 
signals. The noise is introduced by the disturbances and the fluctuations 
beyond the control of the designer. The problem of noise is a problem 
of much research in connection with communications engineering. There 
the cóntrál question is how to devista system efeede so that the effects 
of the unavoidable noise can be minimized and the useful information of 
| diseuse this particular problem noise tillering in € Ë 
the signal is not desta Toyo e problem o is chapler is, 
however, somewhat different: In our problem, the random output is the 
only output of the system. Our purpose in the design of the system, 
particularly the design of the feedback servomechanism. is to obtain with 
a given input an output of the desired statistical characteristics. We shall 
see that the method of transfer function developed in the previous chapters 
remain useful in the present task. 


DI 
X. Statistical Description of a Random Function 


' Let us consider a system which generates a random function ya. 
Y 
Now to formulate the concept of a statistical description of such a random 
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Chapter XII 


Linear System with Variable Coefficients 


The only system with time varying coefficients considered in 
detail in the previous chapters is the pendulum with a periodic force at 
the supporting end, discussed in connection with the phenomenon of 
parametric excitation and damping. All other systems considered do 
not have coefficients of their differential equations that are explicitly 
functions of time. We have, however, shown in Chapter 1, that linear 
systems with time varying coefficients can have behavior entirely different 
from systems with constant coefficients. In this chapter, we shall again 
take up this question and discuss in some detail such a typical but simple 
system: the short range artillery rocket. We shall demonstrate that the 
question of stability of such a system with variable coefficients cannot be 
solved in the same manner as the linear system with constant coefficients. 
Not only is the method of Laplace transform and transfer function useless 
for the purpose, but also we are forced to change our entire approach to 
the problem. 

We shall study the motion of a fin-stabilized artillery rocket 
during the period of action of the rocket thrust. We shall be particularly 
concerned about the angular deviations of the rocket axis from the launch- 
ing angle due to the action of disturbances when leaving the launcher and 
the subsequent damping action of the fins. The general problem of 
dynamics of artillery rockets has been studied in great detail by various 
authors in different countries during World War H. The American work 
is summarized by Rosser, Newton, and Grogs. * The work done in 
England is reported by Rankin, **  Carriere's paper*** represents 
French investigation on the same subject. Our discussion here 
will be greatly simplified and has the purpose of only bringing out 


* J. B. Rosser, R. R. Newton, and G. L. Gross, "Mathematical 
Theory of Rocket Flight", McGraw Hill, New York (1947). 


* R. A. Rankin, Philosphical Transcations, Roy. Soc. of London, (A), 
241:457-585, (1949). 


*** P. Carriére, Mémorial de 1 Artillerie Frangaise, 25: 253-360 (1951). 
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automatic sensing and measuring control system, i.e., an optimalizing 
system which automatically holds the airplane at the measured optimum 
operating conditions. 

Of course, a skilled human operator controls the performance of 
a machine on the optimalizing principle: He watches the instrument 
readings of the inputs and outputs of the machine, and then uses his 
knowledge and experience to decide in what directions should the controls 
be adjusted. The adjusted inputs bring new output readings which have 
to be interpreted by the operator to determine whether the optimum 
operating condition is reached or exceeded. New adjustments of the 
control will have to be made. The continuous adjustment of inputs is 
the sensing process and the reading of the outputs is the feedback. How- 
ever, manually-controlled optimalizing systems are necessarily slow in 
response, and for complicated systems human skill, no matter how 
developed, is not sufficient. Automatic optimalizing control was con- 
ceived by C. S. Draper, Y. T. Liand H. Laning, Jr.* Its application 
to cruise control of airplane was discussed by J. R. Shull. * 


TP d Principles of Optimalizing Control 


The heart of an optimalizing control system is the non-linear . 
component which characterizes the optimum operating conditions. For 
simplicity of discuseion, we shall assume that this basic component has 
a single input and a single output. For the time being, we shall neglect 
also any time effects and assume that the output is determined only by 几 s- nae 


valu 4 the input. Since there is an optimum operating point, output as a function 


of input has a maximum at EI and N as shown in Fig. 15.1. It is 
convenient to refer the output and the input to the optimum point and put 
the physical input as Z WR and the eee output as yr A . The 
optimum point is then the point x m 3 =O. The purpose of an optimal- 


izing control is then to search out this optimum point and to keep the 


* Y. T. Li, Instruments, 25: 72-77, 190-193, 228, 324-327, 350-352 
(1952). C. S. Draper and Y. T. Li, "Principles of Optimalizing 


Control Systems and an Application to Internal Combustion Engine", 
ASME Publications (1951). 1 


** J. R. Shull, Trans. I. R. E. (Electronic Computers), Dec. 1952, 
pp. 47-51. 4 
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we 

optimum filter design be abanioníeg the inadequate RC-circuits and 

by actually using an analog computer or even a digital computer to 

serve as the filter. Then the theoretical optimum performance can 

be actually attained. However, the introduction of electro mechanical 
computer as a component of the filtering system certainly greatly 
increases the complexity of the overall system, and can be justified 
only in very critical cases. But if we have made the system very 
complicated at tire cost, we may ask whether we will actually obtain 

the very best performance. The optimum performance in the theory 
discussed in the previous sections is only optimum within the limitations 
of the assumptions of the theory. For instance, two random signals 
with the same correlation function or the same power spectrum, 
according to the theory developed, require the same optimum filter. 
This is, ina sense, a certain looseness in design criteria. Surely, 

if we have more statistical information about the signal than just power 
spectrum, we should be able to distinguish these two signals and to 
improve our design by utilizing such additional knowledge. Then we 

can obtain even better performance than possible with the so called 
"optimum" filter. It is evident that this generalized approach to the 
filtering problem must require more advanced theory of probability 

than we have used. The recently developed science of information ES A 
theory may also find important applications here. A beginning* has been 
made inthis "probabilistic" approach to problem of detecting signal in 
noise. But much remains to be done. 

The Wiener-Kolmogoroff theory of optimum filter is based upon 
the mean square error criterion. By using this criterion, we essentially 
put the emphasis on minimizing the large errors without much consider- 
ation on the. small errors. However in many occasions, we may be most 
interested in making the frequent error as small as possible, while not 
particular about making an infrequent large error. It is also possible 
that the probability function is very lopsided, with the mean far from 
the mode. For such cases, the mean square error criterion is entirely 


* See for instance, P. M. Woodward, I. L. Davies, Phil. Mag. 
41:1001-1017 (1950); Proc. I. R. E. 39:1521-1524 (1951); J. I. E. E. 
(London) 99(111); 37-51 (1952). T. G. Slattery, Proc, I. R. E. 

40: 1232-1236 (1952). 
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Chapter XVII 


Control of Error 


In the preceding chapter, we have shown how the principle of 
ultrastability can make the control system insensitive to accidental 
errors and occasional failures of the components by the simple device 
of changing the characteristics of the system whenever instability occurs. 
Since an ultrastable system will automatically seek stability, the control 
system, when designed, actually embodies unstable fields of behavior as. 
well as stable fields of behavior. In other words, during the design of 
an ultrastable system, we make no attempt to predict stability from 
instability, to separate the right fields of behavior from the wrong fields 
of behavior. Errors of behavior are merely treated as a probability, 
but otherwise unspecified. In this chapter, we shall approach the 
reliability of complex control system from a different point of view: We 
sball specifically introduce errors into the system and ask how should 
the system be designed so that the system will give satisfactory perfor- 
| mance in spite of the errors. That is, we wish to know how to control 
the error. ; 

This subject of control of error is now in its early period of 
development. Only the control of error in most elementary operation 
can be discussed, And this is wholly due to J. von Neumann. * Our 
discussion in this chapter is then an exposition of Neumann's work. Its 
purpose is to serve as an introduction to this very important topic and 
to indicate the need for much further investigations. 


1f.] Reliability by Duplication 

It is common knowledge that reliability of a system can generally 
be increased by the simple expedient of duplication. For instance, ifa 
simple system as shown in Fig. 18. la has the characteristic that when 
it fails to operate, it merely gives no output. Then to guard against 
probability of failure, we can duplicate the system with % identical 


* J. von Neumann, "Probabilistic Logics and the Synthesis of Reliable 
Organisms from Unreliable Components". Printed notes of lectures É que 
given at the California Institute of Technolo Pasadena, California, 2 
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Egs. (18.27) and (18. 28) have a first term identical with Eq. (18.1). 
The additional terms come from the imperfect elements and from the 
statistical distribution of errors. 

` With any specified E, 7 , € and 2 . Eqs. (18.26), (18.27), 
and (18. 28) enable us to compute the distribution function of J , the 
fraction of activated output lines of the complete Scheffer stroke system. 
We can make this somewhat clearer by reverting to the notation of 
probability distribution functions. Thus for instance Eq. (18, 26) is 
equivalent to 2 
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The probability distribution function of Y , W(7;5,7;2) , is thus the 
result of integrating with respect to g and A of the joint probability 
of 5 , M, and ). Thus 
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We shall now show that under proper conditions we can obtain 
almost perfect performance of the multiplexed Scheffer stroke system 
by increasing 2. Consider a given fiduciary level &. Perfect per- 
formance means the implication of 75 , or non-activation of out- 
put, by E > =í , Giz -, or the activation of both inputs; the 
implication of J2/-S , by either Pai , Y2/-S or F -&, 
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